METHOD OF SPECTRUM ANALYSIS IN TWO-DIMENSIONAL REPRESENTATION 

Field of the Invention 

The present invention relates to a method of spectrum 
analysis in two-dimensional representation to obtain accurate 
characteristic information on an object to be analyzed by using 
derivatives on spectral data of the said object in spectroscopic 
analysis . 

Description of the Prior Art 

Heretofore, in spectroscopic analysis (including 
infrared spectrophotometry and the like) in general, by 
measuring an intensity of light transmitting through the object 
to be analyzed (may hereafter be denoted as "analyzed object") 
and converting absorbance, various characteristics or 
information in terms of physical and/or chemical properties (may 
hereafter be denoted as "characteristic information") including 
shapes, concentration, and the like of substances contained in 
the said analyzed object have been obtained based on the shapes 
of the said spectrum or a spectral profile. As such a 
conventional spectral profile, that spectral profile 
representing absorbance of the analyzed object in terms of the 
wavenumber (or wavelength) are generally used (may hereafter 
be referred to as "absorbance/wavenumber spectral profile") . 
It can be considered that such a spectral profile is formed with 
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a number of component bands overlapped. 

As a means to analyze such spectral profiles, for example, 
the Japanese patent publication Hll-148865 discloses the method 
of using derivatives on the spectral waveform of the analyzed 
object in spectroscopic analysis. To be concrete, as for the 
specific analyzed object, when spectral waveform where the 
intensity of the signal output of a spectrophotometer is 
represented as a function of wavenumber, wavelength, or time 
is prepared letting nl and n2 (nl*n2) be a positive integer, 
a method of spectrum analysis in two-dimensional representation 
with the following characteristics is disclosed. These 
characteristics are; 

Calculating the nl-th and n2-th derivatives with respect 
to wavenumber, wavelength, or time of the said intensity of the 
signal output regarding the spectral profile on the said 
spectral data, plotting points on the two-dimensional 
coordinate plane whose X-coordinate is the said nl-th derivative 
and whose Y-coordinate is the said n2-th derivative, 
respectively, on the said two-dimensional coordinate plane, 
creating a two-dimensional plot of derivative pairs on the said 
spectral waveform, thereby obtaining specific information on 
the said spectral profile based on the two-dimensional plot of 
the said derivative pairs. 

However, although the spectral analysis method cited in 
the Japanese patent publication Hll-148865 can find component 
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bands when the spectral profile contains a (single) component 
band, it is difficult to find all these component bands and even 
to estimate them is extremely difficult when the spectral profile 
of the analyzed object contains several component bands 
overlapped. In particular, there are many cases that the 
spectral profile of the analyzed object has several overlapping 
component bands. Therefore, as for the analyzed object with a 
spectral profile containing several component bands, the method 
described in the Japanese patent publication Hll-148865 can 
estimate the dominant component bands to some extent among the 
said component bands. However, as for unresolved bands buried 
under the dominant band, that is, the characteristics of those 
component bands are unclear due to severe overlap, it is 
extremely difficult even to estimate those bands. 

SUMMARY OF THE INVENTION 
Based on the method cited in the Japanese patent 
publication Hll-148865, the present invention has been 
developed for further improvement. As a result of intensive 
studies, as for the specific analyzed object, when spectral data 
where the intensity of the signal output of a spectrophotometer 
is represented as a function of wavenumber, wavelength, or time, 
are prepared letting n and m (n^m) be a positive integer, the 
inventor has proposed a method of spectrum analysis in two- 
dimensional representation for obtaining the specific 



3 



characteristic information on the said spectral data based on 
the said two-dimensional derivative plot by calculating the n-th 
and m-th derivatives with respect to wavenumber, wavelength, 
or time, of the said spectral data, plotting points on the 
two-dimensional coordinate plane as the X-Y coordinate system 
whose X-coordinate is the said n-th derivative and whose Y- 
coordinate is the said m-th derivative respectively on the said 
two-dimensional coordinate plane, and preparing the two- 
dimensional derivative plot on the said spectral data. 

Further, based on the characteristic information 
described above, the inventor has proposed a method of spectrum 
analysis in two-dimensional representation which estimates the 
component bands comprising the spectral profile of the analyzed 
object by estimating band parameter values regarding at least 
one component band among the component bands contained in the 
spectral profile of the analyzed object, estimating at least 
one component band, obtaining the two-dimensional derivative 
plot with the specific remaining component bands removed by 
clearing the specific component band or specific component bands 
already estimated or the two-dimensional derivative plot from 
spectral profiles or the two-dimensional derivative plot of 
analyzed object, obtaining specific characteristic information 
based on the two-dimensional derivative plot of this specific 
component removed, estimating band parameter values on other 
component bands based on the said characteristic information, 
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and iterating the estimation of at least one of the other 
component bands thereby estimating component bands in order. 

That is, the present invention is a method of estimating 
component bands in order, using a specific component band or 
specific component bands already estimated and estimating 
remaining component bands (including the ones not yet estimated, 
the ones desired to be optimal) . 

To be concrete, a preferred embodiment of the present 
invention, although not specifically limited, is, first of all, 
to estimate one component band BDi(i is a positive integer) by 
estimating the band parameter values for BDi among component 
bands contained in the spectral profile of the analyzed object 
based on the specific characteristic information on the spectral 
profile of the two-dimensional derivative plot. 

Next, the preferred embodiment is to prepare the two- 
dimensional derivative plot with BDi removed, either creating 
the two-dimensional derivative plot with BDi removed by clearing 
the two-dimensional derivative plot of the component band BDj ( j 
is a positive integer where j^i) from the two-dimensional 
derivative plot of the said analyzed object or obtaining the 
spectral profile with BDi removed by clearing a profile of BDi 
from a spectral profile of the said analyzed object. 

Then, the preferred embodiment is to estimate one 
component band BDj by estimating band parameter values on BDj 
other than BDi among the component bands contained in the 
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spectral profile of the analyzed object based on the said 
characteristic information by obtaining the specific 
characteristic information based on the said two-dimensional 
derivative plot with BDi removed. 

Next, the preferred embodiment is to prepare a two- 
dimensional derivative plot with BDj removed, either creating, 
the two-dimensional derivative plot with at least BDj removed 
by clearing BDj or the two-dimensional derivative plot of both 
BDi and BDj from the two-dimensional derivative plot of the said 
analyzed object or obtaining the spectral profile with at least 
BDj removed by clearing BDj or the profile of both BDi and BDj 
from a spectral profile of the said analyzed object. 

Then, the preferred embodiment is to estimate the 
component band BDk (k is a positive integer where at least k* j ) 
by estimating band parameter values on BDk other than BDj among 
the component bands contained in the spectral profile of the 
analyzed object based on the said characteristic information 
by obtaining the specific characteristic information based on 
the said two-dimensional derivative plot with BDi removed. 

Thus, the present invention is a procedure of estimating 
other component bands in order as for component band contained 
in the spectral profile of the analyzed object, using the 
specific component band or specific component bands already 
estimated by iterating estimation of component bands, creating 
the two-dimensional derivative plot with information on an 
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estimated component band or estimated component band removed 
from the information on the analyzed object, acquiring the 
specific information based on this two-dimensional derivative 
plot and operating the estimation of other component bands based 
on this characteristic information. 

In the two-dimensional derivative plot on spectral data 
of the analyzed object, the characteristic information on 
dominant component bands can be distinguished to some extent, 
however, as for characteristic information on unresolved 
component bands buried under the dominant component band, since 
it is subtle or since it overlaps with the characteristic 
information of other component bands overlapping with the said 
component band, the characteristic information of the said 
component band of its own becomes unclear and it is difficult 
to be distinguished. 

Particularly, in the case where the overlapping component 
bands are dominant ones, the characteristic information on the 
unresolved component bands buried under the dominant component 
band can seldom be distinguished. 

However, as the analytic method of the present invention 
shows, when the two-dimensional derivative plot is created by 
removing information on dominant component bands estimated from 
the characteristic information of the two-dimensional 
derivative plot of the analyzed object from the information about 
the analyzed object and by obtaining the two-dimensional 
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derivative plot based on the spectral profile after removal of 
information on dominant component bands, the characteristic 
information of component bands hidden by overlapping dominant 
component bands appears, which enables to obtain the 
characteristic information about component bands hidden by 
overlapping, and further enables to obtain the characteristic 
information which has never been achieved by the two-dimensional 
derivative plot on the spectral data of the analyzed object. 

Then, by iterating the operations of removing the profile 
of single or several component bands already estimated from the 
information about the analyzed object and by obtaining the 
characteristic information using the two-dimensional 
derivative plot based on the removal thereafter estimating other 
component bands, estimating the component band contained in the 
spectral profile of the analyzed object in order with suitability 
can be achieved. 

Thus, the method of estimating a single component band 
or several component bands contained in the spectral profile 
of the analyzed object and finding the other component bands 
in order by using the two-dimensional derivative plot obtained 
by removing the profile of the single component band or several 
component bands from the spectral profile of the analyzed object, 
or this procedure is called "Band Stripping". 

Therefore, when the method of the present invention is 
employed, the component band that is contained in the spectral 
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profile of the analyzed object can be estimated easily. 

Further, in the preferred embodiment of the present 
invention, the component band is a Gaussian band, a Lorentzian 
band, or the mixture thereof. 

In the present invention, n is preferably 1 and/or 3 and 
m is n + 1. When the combination of n and m is (n,m) = (1,2), (3, 4) , 
the characteristic information on the spectral data appears 
clearly. Therefore, when the two-dimensional derivative plot 
is created with the combination of (n, m) = ( 1 , 2 ) and/or (3,4), the 
characteristic information can be obtained easily. 

According to the method of spectrum analysis in two- 
dimensional representation, in the two-dimensional derivative 
plot where pairs of the first and second derivatives are 
represented in the X-Y coordinate system, when a typical local 
minimum indicates the existence of a corresponding component 
band, an X position of the said local minimum is a first 
approximation of band center position Xc of the said component 
band, setting several points on the said two-dimensional 
derivative plot in the vicinity of Pd, point of intersection of 
the said two-dimenstional derivative plot with the X- axis, as 
candidates for the inflection point of the said component band, 
estimating the bandwidth of the said component band from the 
candidate of the said inflection point by the following Equation 
(1), estimating the peak height of the said component band from 
the distances between the said local minimum and the point (s) 
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of intersection of vertical line passing through the said local 
minimum and the horizontal line(s) passing through the said 
candidate points, obtaining the candidates for band parameter 
values of the said component band, and further obtaining the 
constraint conditions subjected to the band parameter values 
for the said component band from the said two dimensional 
derivative plot, the relation between the bandwidth bw and the 
X-position of the inflection point X P of a single band can be 
preferably expressed by 

bw = (1/Kp) lXc- XpI (1) 
(In Equation, bw is an estimated value of the bandwidth 
of a Gaussian or a Lorentzian band, where the coefficient Kp is 
0.42466 for Gaussian and 0.288675 for Lorentzian.) 

According to the method of spectrum analysis in two- 
dimensional representation, in the two-dimensional derivative 
plot where pairs of the third and fourth derivatives are 
represented in the X-Y coordinate system, when a typical local 
maximum indicates the existence of a corresponding component 
band, an X position of the said local maximim is a first 
approximation of band center position Xc of the said component 
band, setting several points on the said two-dimensional 
derivative plot in the vicinity of Qd, point of intersection of 
the said two-dimenstional derivative plot with the X- axis, as 
candidates for the secondary inflection point of the said 
component band, estimating the bandwidth of the said component 
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band from the candidate of the said secondary inflection point 
by the following Equation (2) , estimating the peak height of 
the said component band from the distances between the said local 
maximim and the point (s) of intersection of vertical line passing 
through the said local maximum and the horizontal line (s) passing 
through the said candidate points, obtaining the candidates for 
band parameter values of the said component band, and further 
obtaining the constraint conditions subjected to the band 
parameter values for the said component band from the said two 
dimensional derivative plot, the relation between the bandwidth 
bw and the X-position of the secondary inflection point Xq of 
a single band can be preferably expressed by 

bw = (1/Kp) |Xc- Xq| (2) 
(In the Equation, bw is an estimated value of the bandwidth 
of a Gaussian or a Lorentzian band, where the coefficient Kq is 
0.31508 for Gaussian and 0.16426 for Lorentzian.) 

In the present invention, the method of spectral analysis 
in two-dimensional representation which adjusts the already 
estimated band parameter values can preferably be used so that 
the specific component band already estimated and the 
complementary estimated component band with all the estimated 
component bands other than the said estimated specific component 
band removed from the spectral profile of the analyzed object 
or two-dimensional derivative plot coincide. 

With the said band stripping method, the parameter values 
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of the each component band cannot be truly estimated and it 
sometimes fails into a so-called local optimum. This is because 
the adjacent component bands overlap and when the parameter 
values of the specific component band are too large or too small, 
the overestimation and/or underestimation affect (s) the 
parameter values such as band center position, bandwidth, and 
peak height of the adjacent component bands. Therefore, the 
present invention is capable of estimating the band parameter 
values for the component band with further suitability by 
introducing "complementary estimated component band". 

The example of the coincidence between the specific 
estimated component band (denoted as eBD) and its complementary 
estimated component band (cBD) includes a method of minimizing 
the total sum of the distance of the iso-wavenumber lines between 
eBD and cBD in the two-dimensional derivative plots. By 
adjusting the parameter values for the specific estimated 
component band and by minimizing the difference between eBD and 
cBD, the estimated parameter values become as close to the true 
value as possible. In addition, by improving the degree of 
symmetry of the complementary estimated component band, the 
adjacent band parameter values are optimized. In other words, 
by adjusting the band parameter values for the specific estimated 
component band, the degree of symmetry of the complementary 
estimated component band is improved. Thus, introducing the 
complementary estimated component band and aiming to optimize 
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the estimated band parameter values already obtained is called 
"Complementary Matching" method. 

The object of the present invention is to provide a method 
of spectrum analysis that can easily estimate its several 
component bands of the analyzed object having spectral profile 
containing several component bands. In addition, it can be 
applicable to spectral data such as infrared spectra, visible 
light spectra, ultraviolet spectra, Raman spectra, X-ray 
dif f ractogram, and chromatogram, and the like. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a graph depicting an original spectrum having 
three component bands overlapped. 

Figure 2 is a graph depicting a second derivative spectrum 
of the original one shown in Fig.l. 

Figure 3 is a graph depicting a fourth derivative spectrum 
of the original one shown in Fig.l. 

Figure 4 is D1-D2 plots for 0V3 and OV12 . 

Figure 5 is D3-D4 plots for OV3 and OV12 . 

Figure 6 is a graph showing the permissible range of the 
bandwidth and peak height for BD3 estimation. 

Figure 7 is D1-D2 plots for OV12 and BDi . 

Figure 8 is D1-D2 plots for OV13 and BDi. 

Figure 9 is D1-D2 plots for OV13 and BD3. 

Figure 10 is D1-D2 plots for OV12 and eOVi2 . 

Figure 11 is a graph showing the permissible range of the 
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bandwidth and peak height for BD2 estimation. 
Figure 12 is D1-D2 plots for eOVi3. 

Figure 13 is a graph showing the permissible range of the 
bandwidth and peak height for BD3 estimation. 

Figure 14 is D1-D2 plots for BDi and eBDi . 

Figure 15 is a graph depicting spectral profile of 
eBDi, eBD2, eBD3, eOV3, and OV3. 

Figure 16 is D1-D2 plots for eOV3 and OV3. 

Figure 17 is D1-D2 plots for eBD3 and cBD3 . 

Figure 18 is D1-D2 plots for eBD 2 and CBD2. 

Figure 19 is D1-D2 plots for eBD 1 and eBDi. 

Figure 20 is D1-D2 plots for eBD 3 and cBD3 after the 
parameter values are improved by SumLS . 

Figure 21 is D1-D2 plots for eBD 2 and CBD2 after the 
parameter values are improved by SumLS. 

Figure 22 is D1-D2 plots for eBD 1 and eBDi after the 
parameter values are improved by SumLS. 

Figure 23 is D1-D2 plots for eBD 3 and cBDs after the 
parameter values are improved in coincidence and symmetry. 

Figure 24 is D1-D2 plots for eBD 2 and CBD2 after the 
parameter values are improved in coincidence and symmetry. 

Figure 25 is D1-D2 plots for eBDi and eBDi after the 
parameter values are improved in coincidence and symmetry. 
DESCRIPTION OF THE PREFERRED EMBODIMENT 

The present invention is explained in detail based on the 
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drawings. Figure 1 shows the original spectrum generated (or 
synthesized) that is consisted of the three overlapping 
component bands. The explanation goes as for a method of 
estimating component band from these graphs. 

As a shape of the component band estimated by the spectral 
analysis method of the present invention, a Gaussian band, a 
Lorentzian band, or a mixture thereof can be considered. 

A Gaussian band is expressed by the following Equation: 
BDg(X)= phc exp {-4 log 2 (X-bcc) 2 /bwc 2 } (3) 
(in the Equation, bcc is the center position of the band, phc the 
peak height, and bwc2 the bandwidth.) 

A Lorentzian band is expressed by the following Equation: 
BDl(X)= ph L / {1 + 4 (X-bCL) 2 /bwL 2 } (4) 
(in the Equation, bcL is the center position of the band, phL the 
peak height, and bwL the bandwidth.) 

A mixture of a Gaussian and a Lorentzian band is expressed 
by the following Equation: 

BDm(X)= mBDc(X) + (l-m)BDL(X) (5) 
(in the Equation, BDc(X)is the Gaussian band given by Equation 
(3) , BDl(X) the Lorentzian band given by Equation (4) , and m is 
the ratio of the mixture where 0<m<l.) 

Further, the bandwidth is the full-width at half-height 
of the band (FWHH) . In other words, the bandwidth is the distance 
between the two X positions at half-height of the band. 

In the present invention, in order to find one component 
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band which is contained in the original spectrum, the 
characteristic information is so obtained as by differentiating 
the original spectrum. M n-th-m-th two-dimensional derivative 
spectral profile'' is so defined as the one with the n-th 
derivative on the X-axis and the m-th derivative on the Y-axis. 
In some cases, it is abbreviated as xx Dn-Dm plot". 

Figure 1 shows the original spectrum of 0V3 (=BDi+BD2+BD3) 
having BDi, BD2, and BD3, three Gaussian component bands 
overlapped and spectral profile of each component band. That 
is, OV3 is obtained by summing the component bands calculated 
by the following Equation (6) . 

BDi(X)=phi exp {-4 log 2 (X-bci) 2 /bwi 2 } (6) 
(in the Equation, i is 1, 2 or 3, representing the i- 
th band, bci the center position of the band, phi the peak height, 
and bwi the bandwidth.) 

Spectral analysis is performed with Mathematical 
mathematical software packages (ver. 2. 2, Wolf ram Research, Inc., 
IL) . The band parameter values for the individual component 
band are given as follows: 

Component band BDi : bci = -0.5978 

bwi = 1.1070 
phi = 0.6504 
Component band BD2: bc2 = 0.0012 

Bw2 = 0. 9643 
ph2 = 1.0630 
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Component band BD3 : bc3 = 0.4533 

Bv/3 = 0.8078 
pha = 0.8123 

In the simulation study estimating the component band, 
the parameter values are treated as unknown, however, the true 
parameter values described above are sometimes utilized in the 
course of the following explanation for evaluation of obtained 
results . 

Mathematica is also used for calculations and drawings. 
In differential operation, analytic differentiation is 
performed and digitized with an interval of 0.01 so as to avoid 
the error caused by digital differentiation. 

Here, in the usual spectral profile, the X-coordinate 
represents wavenumber, wavelength, or time, whereas the Y- 
coordinate represents absorbance. For convenience, the 
explanation goes defining that the X-coordinate represents 
wavenumber and the Y-coordinate represents absorbance. Further, 
the range of the X-coordinate is -2 < x < 2. Here, in this 
specification, "the n-th derivative spectrum" is denoted as the 
derivative spectrum with the guantity of the Y-axis 
differentiated the n-th with respect to the X-axis . In addition, 
"two-dimensional derivative plot" is denoted as the derivative 
plot with the n-th derivative on at least one axis of the two 
axes in the two-dimensional coordinate plane as the X-Y 
coordinate system. 
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In the analytical procedure of the method of this invention, 
first, a two-dimensional plot of the analyzing spectrum is 
created, and band parameter values for at least one component 
band within the two-dimensional plot are estimated after 
obtaining the characteristic information. For this purpose, 
the characteristics of the component band are explained as below. 
(Characteristic points of the component band) 

In Fig. 1, among the three component bands, the 
characteristic points are shown as for the sharpest (with the 
narrowest bandwidth) component band BD3. Point T3 (Xt3=0.4 533) 
is the vertex. (Xt3 represents the X position of the point T3 
and the same representation follows hereafter.) Points P31 
(Xp3i=0 . 1103) and P32 (Xp32=0 . 7 963 ) are the zero-crossing points 
of the second derivative and the inflection points. Points Q32 
(Xq32=0. 1988) and Q33 (Xq33=0 . 707 8 ) are the innermost zero- 
crossing points of the fourth derivative and the inflection 
points. Hereafter, the innermost zero-crossing points of the 
fourth derivative are called secondary inflection points. 

In estimating the parameter values for component band BD3, 
band center position bc3 can be estimated by the corresponding 
local minimum or maximum, and bandwidth bw3 can be estimated by 
the following Equation. However, peak height ph3 cannot be 
estimated easily. When point Bk is so placed as the point where 
X=Xt3 on OV12 which is obtained by removing BD3 from OV3, and when 
point TV is so placed as the point where X= Xt3 on OV3, the length 
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of line segment T3'Bk corresponds to peak height ph3. Therefore, 
if the X position of point Bk can be estimated, peak height ph3 
can also be estimated. 

Bandwidth bwi of component band BDi can be calculated from 
Xpii, the X position at inflection point Pii, or from Xp±2, the 
X position at inflection point Pi2, by the following Equations: 
bwi = (1/Kp) |bci-xpii| (7) 
bwi = (1/Kp) |bci-x P x2 | (8) 
(Kp= 0.4247 in the case of a Gaussian band, whereas Kp = 0.2887 
in the case of a Lorenztian band.) 

Bandwidth bwi of component band BDi can be calculated from 
Xq±2, the X position at the secondary inflection point Qi2, or 
from Xqi3, the X position at secondary inflection point Q±3, by 
the following Equations: 

bwi - (1/Kq) Ibci-xoii | (9) 
bwi = (1/Kq) |bci-XQi 2 | (10) 
(Kq= 0.3151 in the case of a Gaussian band, whereas Kq = 0.1625 
in the case of a Lorenztian band.) 

In order to estimate each component band, it is necessary 
to find a typical local minimum, local maximum, and the 
inflection points of the corresponding component band on OV3 
where all component bands overlap. For this purpose, first, the 
feature of OV3 is studied by means of its derivative. 

The second derivative spectrum of the original one in Fig. 
1 is shown in Fig. 2, while the fourth derivative spectrum in 
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Fig.l is shown in Fig. 3. OV12 is the spectrum where BD3 is 
subtracted from OV3 (i.e. , OV12 = 0V3 - BD3 ) or the spectrum where 
BDi and BD2 overlap ( i . e . , OV12 = BDi + BD2 ). In Fig. 2, Min3 is 
the local minimum, points Pa and Pd are the inflection points 

of 0V3. 

Points P31' and P32' correspond to inflection points P31 
and P32 of BD3 and they are the isosbestic points in the second 
derivative. As in the initial spectral profile, the length of 
line segment T3'Bk corresponds to the peak height of BD3. It 
should be noted that a peak height of a differential spectrum 
does not equal that of an original spectrum, and it should be 
calculated based on a peak height at the second-order 
differentiation of a Gaussian band with a unit peak height (phi=l) 
having the same bandwidth. 

In the fourth derivative in Fig. 3, Max3 is the local 
maximum, and points Qc and Qd are the inflection points of the 
second derivative of OV3 . Points Q32' and Q33' correspond to 
secondary inflection points Q32 and Q33 of BD3 and they are the 
isosbestic points in the fourth derivative. The length of 
segment TVBk corresponds to the peak height of BD3. 

It should be noted that Min3 is the local minimum of the 
second derivative spectrum of OV3, and Max3 is the local maximum 
of the fourth derivative spectrum of OV3. A typical local 
minimum or local maximum indicates that a corresponding 
component band exists in the vicinity thereof. According to the 
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illustration of Fig. 2, only Min3 is the local minimum of the 
second derivative spectrum of OV3 and there is no local minimum 
that corresponds to other component bands BDi and BD2. On the 
other hand, according to the illustration of Fig. 3, the typical 
Max3 exists as the local maximum of the fourth derivative 
spectrum of OV3. Further, although there are other two local 
maxima, since they are modulated by the side lobe of the dominant 
peak, it is not clear whether they indicate the existence of 
the corresponding component bands. 

From the description mentioned above, BD3 is the analyzing 
component in which the local minimum is typically at the second 
derivative spectrum and in which the local maximum is typically 
at the fourth derivative and the remaining components OV12 is 
called "background components''. Then, estimating band 
parameter values is examined based on the algebraic geometry. 
In a spectrum where many component bands overlap, the 
generalization is possible by regarding it as the two-component 
system of the "analyzing" component , band and the "background" 
components. "Two-dimensional derivative plot" refers to 
plotting the pairs of the n-th and m-th derivatives in the 
two-dimensional coordinate. It can be abbreviated as "Dn-Dm" 
plots . 

D1-D2 plots are shown in Fig. 4, in the two-dimensional 
representation of the pairs of the first derivative and the 
second derivative. D3-D4 plots are shown in Fig. 5, in two- 
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dimensional representation of the pairs of the third derivative 
and the fourth derivative. 

Letters and sequential numbers are attached as for the 
characteristic points of the corresponding component band BD3. 

That is, in Figs. 4 and 5, points T3' , P31' , P32' , Q32' , and 
Q33' on OV3 correspond to points T3, P3i,P32,Q32, and Q33 on BD3, 
respectively. 

In D1-D2 plots in Fig. 4, if the position of T3' , P31' , or 
P32' can be estimated, or in D3-D4 plots in Fig. 5, if the position 
of T3' , Q32' , or Q33' can be estimated, the band center position 
of BD3 and its bandwidth can be determined by Equation (3) or 
(4) . 

In D1-D2 plots in Fig. 4, the local minimum Min3 can be 
a candidate for T3' , and the point of intersection of Dl-axis, 
Pd, can be a candidate for P32' . However, it should be noted that 
the X position of Min3 (Xmin3=0 . 3285 ) is away from that of point 
T3' (Xt3=0 . 4533) , the true value, and the X position of point Pd 
(Xpd=0.6746) is far away from that of point P32' (Xp32=0 . 7963) . 

On the other hand, in D3-D4 plots in Fig. 5, the local 
maximum Max3 can be a candidate for TV , and Qd, the point of 
intersection of D3-axis can be a candidate for Q33' . X position 
of Max3 (Xma*3=0 . 4180) is closer to that of point T3' . Further, 
X position of point Qd (Xod=0.6482) is closer to that of point 
Q33' (Xq33=0.7078) , the true value. Concerning the two- 
dimensional derivative plot, since D1-D2 plots have simpler 
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shapes than do D3-D4 plots, the explanation goes transferring 
the characteristic points on D3-D4 plots to those of D1-D2 plots, 
where the point corresponding to Max3 is denoted Max3' and the 
inflection point converted by the secondary inflection point 
Qd is denoted Pa' . 

Iso-wavenumber lines are used for explaining the 
geometrical relation of the characteristic points on the 
two-dimensional derivative plot, Iso-wavenumber lines are the 
straight lines connecting the same wavenumber points between 
two spectra. In Figs. 4 and 5, since it is complicated to draw 
the lines over entire range, the iso-wavenumber lines are drawn 
only in the region of interest. Here, the sequence numbers along 
the profile represent the wavenumbers . 

In Fig. 4, the iso-wavenumber lines passing through point 
T3' is parallel to D2-axis. Let the point of intersection with 
OV12 be point Bk, let the point of intersection of the extended 
iso-wavenumber line and OV3 be point At. The length of line 
segment Ta'Bk corresponds to the peak height of BD3. 

On the other hand, the iso-wavenumber lines passing 
through the inflection points P32' and P31' , respectively, are 
parallel to Dl-axis. Let the point of intersection of the 
iso-wavenumber lines passing through point P32' and OV12 be point 
Pj, and let the point of intersection of the extended iso- 
wavenumber line with respect to P32' and OV12 be point A P . Let 
the point of intersection of the iso-wavenumber line passing 
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through point P31' and OV12 be point Pi . The length of line segment 
P32'Pj equals the length of line segment P3i'Pi, which corresponds 
to a half of the bandwidth between two inflection points of BD3. 
Let the point of intersection of lines T3'Bk and P32'Pj be point 
Bj, and let the point of intersection of lines T3' Bk and P32' Pi 
be point Bi. Further, let the point of intersection of lines 
T3'Bk and P31' P32' be point Bx. 

According to Fig. 4, since points Bk and Bx are above the 
line which is parallel to Dl-axis passing through point Pd and 
are close to each other, it causes no problem when point Bk is 
replaced by point Bx. In general, the position of point Bk depends 
on the shape of background component OV12 and it is difficult 
to find its position by the algebraic geometrical method. 

Nevertheless, the permissible region of point Bk can 
approximately be determined. When the component bands are close 
to each other, the shape of OV3 is like that of a single band. 
At this time, the position of point Bk is below that of point 
Bx, located in the vicinity of point Bi. Conversely, as the 
separation between the component bands is large, the position 
of point Bk moves upward to that of point Bx. Further, when the 
separation between the component band becomes larger, another 
local minimum appears on OV3 and the position of point Bk is 
located in the vicinity of point Bj which is above point Bx. 

On the other hand, as for line segment P32' Pj corresponding 
to a half of the bandwidth between two inflection points of BD3, 
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point Pj exists in the left-hand side of point A P . That is, line 
segment P32'Pj is shorter than line segment P32'A P . Likewise, the 
same analysis can be applied to D3-D4 plot in Fig. 5. Based on 
the geometrical study described above, how to find estimated 
parameter values of a component band is explained one by one. 
(Stepl) 

As can be seen from Fig. 4, the X position of Max3' (or 
Min3) is the first approximation of band center position Xc for 
the corresponding component band BD3. Hereafter this point is 
denoted as eT3. Draw a straight line LI parallel to Dl-axis 
passing through point eT3. An estimated point eP32, on OV3, for 
the inflection point P32 of BD3 is placed in the vicinity of the 
inflection point Pa of OV3. In more detail, the X position of 
eP32 is at the positive side of that of point Pd. Then, an 
estimated value of the bandwidth of BD3 is calculated by Equation 
(8) . Draw a straight line L2 parallel to D2-axis passing through 
point eP32. Line L2 intersects line LI at point Bj . 

Next, an estimated point eP3i, on OV3, for another 
inflection point of BD3 is placed at the opposite side of point 
eP32 with respect to point eT3. Both points eP3i and eP32 are 
equidistant from point eT3 with respect to the X-coordinate or 
the wavenumber. Line segment eP32 eP3i intersects line LI at 
point Bk. From the length of line segment eT3Bj or eT3Bk, an 
estimated value of the peak height of BD3 is calculated. 

Thus, a set of estimated values of band center position 
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eBC3, bandwidth eBW3, and peak height ePH3 for BD3 can be obtained. 
In order to obtain the true values or optimal estimated values 
for BD3, a candidate can be created by finding sets of parameter 
values on several sequential points as the estimated value of 
eP32 from the vicinity of point Pa. 

It is clear that when an estimated value of band center 
position eT3 is close to point T3' , a better estimated value of 
a bandwidth can be obtained. Table 1 shows the sets of estimated 
values of the bandwidth and the peak height for BD3 calculated 
by using systematic sequential points for eP32 when point eT3 
is equal to the true target point T3' . eBW3 is the estimated 
bandwidth, ePH3A is the peak height obtained by line segment 
eT3Bj , and ePH3B is the peak height obtained by line segment 
eT3Bx. 
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Table 1 



Set of estimated value 


X coordinate of eP32 


eBW 3 


ePH3A 


ePH3B 


EstBD3a 


0.70 


0.581 


0.357 


0.244 


EstBD3b 


0.71 


0.604 


0,404 


0.283 


EstBD3c 


0.72 


0.628 


0.455 


0.325 


EstBD3d 


0.73 


0.652 


0.509 


0.372 


EstBD3e 


0.74 


0.675 


0.566 


0.423 


EstBD3f 


0.75 


0.699 


0.628 


0.478 


EstBD3g 


0.76 


0.722 


0.693 


0.537 


EstBD3h 


0.77 


0.746 


0.761 


0.602 


EstBD3i 


0.78 


0.769 


0.833 


0.671 


EstBD3j 


0.79 


0.793 


0.908 


0.744 


EstBD3k 


0.80 


0.816 


0.986 


0.823 


EstBD31. 


0.81 


0.840 


1.07 


0.907 


EstBD3m 


0.82 


0.863 


1.15 


0.995 


EstBD3n 


0.83 


0.887 


1.24 


1.09 
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Figure 6 shows the plots of the candidates of estimated 
values listed in Table 1 with a bandwidth at the Horizontal axis 

and with a peak height at the Vertical axis. Points a, b, c, , 

n on the curve xy (this is called SelLNl) are calculated when 
peak heights are estimated by the length of line segment eT3Bj 

and sequence a', b' , c' , , n' on the curve x'y' (this is called 

SelLN2) are calculated when peak heights are estimated by the 
length of line segment eT3Bx. The point indicated as 'actual' 
is the true point for BD3. It is necessary to restrict estimated 
values from the algebraic geometrical constraint conditions. 

According to Fig. 6, the constraint conditions subjected 
to the estimated values (set of the bandwidth and the peak height) 
are explained as below. As clearly shown in the spectral profile 
of Fig.l, since BD3 must be downward of 0V3, it must be 
BD3 (X) <OV3 (X) . Since BD3 must not exceed OV3 at the estimated 
point of eT3, the band center position of BD3, it must be 
BD3 (eT3) <0V3 (eT3) . There is the upper limit for the peak height. 
That is, in Fig. 6, estimated values must lie below the straight 
line shown as 'OverTop' . 

In addition, estimated values must be at the left-hand 
side of the curved line denoted as 'OverEnv' . When the bandwidth 
is small, a certain degree of the peak height is permissible, 
however, when the bandwidth is large, only a relatively small 
range of the peak height is permissible. Therefore, as shown 
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in Fig. 6, qualitatively, OverEnv is the declined curved line 
from left to right. 

Further, from D1-D2 plots shown in Fig. 4, the constraint 
conditions subjected to the estimated values can be calculated. 
From Fig. 4, the position in which the inflection point P32' of 
BD3 will exist must be above point Pd, the inflection point of 
OV3. The inflection point for an estimated component band is 
above point Pd, which means that the bandwidth of the estimated 
component band is larger than that calculated from point Pd. 
Correspondingly, this means that the estimated value must be 
located at the right-hand side of line PDw in Fig. 6. In addition, 
likewise, since the points Bj and Bx are above point Pd, line 
segment eT3Bj or eT3Bx ought to be longer than the distance of 
eT3 from the X-axis. Correspondingly, therefore, in Fig. 6, 
estimated values must be above line PDh. 

Since point Pd' in Fig. 4 which corresponds to inflection 
point Qd in Fig. 5 has an effect of higher order derivatives, it 
is clear to be generally closer to point P32' . Therefore, 
likewise, in Fig. 6, the estimated values must be located above 
line PDh' and right-hand side of line PDw' . In addition, from 
point At which was already described, constraint of the peak 
height is introduced, which means in Fig. 6 that the estimated 
values are located below straight line OverAt. Therefore, the 
constraint conditions subjected to a bandwidth and a peak height 
or permissible range thereof, are within the region of the 
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hatched polygon bounded by straight lines PDw' , PDh' , OverEnv, 
and OverAt. 

Which point should be selected as the estimated values 
from selection line SelLNl or SelLN2 is explained below. At x' 
end side of SelLN2, both the bandwidth and the peak height are 
smaller than the true value. Conversely, since y' end side of 
SelLN2, both the bandwidth and the peak height are larger than 
the true value and they exceed the curved line OverEnv. Points 
j' , k' , and 1' in the vicinity of the middle part of SelLNl are 
located close to the true value for BD3. The upper limit of the 
bandwidth using point A P (the length of line segment P32' Pj 
corresponding to a half of the bandwidth of BD3 is smaller than 
that of line segment P32'A P ) is in the middle of points k' and 
1' . The position of the true value for BD3 is around the center 
of the hatched triangle region shown in Fig. 6 when the degree 
of overlapping among the component bands is high. And the 
position of the true value for BD3 is around the curved line 
OverEnv when the degree of overlapping is low. Around point k' 
is a good candidate as the estimated values for eBW3 and ePH3. 

The procedure of finding estimated values for BD3 was 
explained as mentioned above, assuming that the estimated value 
of the band center position is the true one. In general, the 
estimated value of band center bc3 of component band BD3 is 
selected as Max3' (Xmax3=0 . 4 18 ) . 

Then, finding a set of a series of estimated values eBD3, 
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estimated values for BD3 are selected. At this stage, a set of 
estimated values for BD3 is obtained (€603=0.418, eBW3=0.830, 
and ePH3=0.780) : 

Thus far, as for BD3, eBW3, and ePH3, estimated values for 
three parameters, eBC3, eBW3, ePH3, were obtained. Here, the 
shape of the component band is assumed to be Gaussian. If the 
shape of the component band is assumed to be Lorentzian, since 
the bandwidth is 1.47 times as large as that of Gaussian, the 
component exceeding OV3 appears. Therefore, the Lorentzian 
component band is rejected. 

After component band BD3 was estimated, a band stripping 
operation is performed, further estimating a component band. 
In order to explain this procedure, as an ideal case, algebraic 
geometrical explanation goes. Figure 7 shows overlaid D1-D2 
plots of OV12 (= BD1+BD2) wherein two component band BDi and BD2 
overlapped and BDi. Figure 8 shows D1-D2 plots of OV13 (= BD1+BD3) 
wherein two component bands BDi and BD3 overlapped and BDi. 
Figure 9 shows D1-D2 plots of OV13 and BD3. 

As can be seen from Fig. 7, since Min2 and Max2' are not 
close to each other, it is more likely or apparent that overlap 
between the component bands BDi and BD2 is strong. 

In Figs. 8 and 9, the local maximum Maxl3 indicates that 
the typical valley exists between BD3 and BDi. The appearance 
of the valley arises from the lower degree of overlap compared 
with that of OV12. In Fig. 8, the curved portion from the origin 
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O to Mini via point Pa substantially overlaps with the curved 
portion of the right half of BDi. Since Mini and Maxl' are not 
close to each other, around the vertex of BDI, they are affected 
by overlap of BD2 and/or BD3. In Fig. 9, the curved portion from 
the origin 0 to Min3 via point Pa substantially superposes on 
the curved portion of the left half of BD3. Since Min3 and Max3' 
are close to each other, the band center position of BD3 is in 
the vicinity thereof. 

In Fig. 8, when estimating for BD3, there. is no problem 
when point Bx is used for estimating point Bk. However, since 
there is the valley shown in Fig. 9, the curved portion of the 
left-hand side of BDi is greatly deformed and since line segment 
Pi2'Pn' and the perpendicular line to Dl-axis passing point TV 
do not intersect, point Bk exists at a position with line segment 
P12' P11' extended. In addition, the curved portion from the 
origin 0 to Min 1 via point Pa substantially superposes on the 
curved portion of the right half of BDi and point Bk is in the 
vicinity of point Bx. Therefore, it is desirable to calculate 
the peak height from the length of line segment Ti'Bi. 
(Step 2) 

Now that estimated values were obtained on BD3, BD3 can 
be expressed mathematically with Equation (6) . This is denoted 
as eBD3. eOVi2 is found by band stripping eBD3 from OV3. That 
is, eOVi2 = OV3-eBD3. In Fig. 10, D1-D2 plots are shown regarding 
OV12 and eOVi2. For comparison, the iso-wavenumber lines are also 
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drawn. Band parameter values for new component band eBD2 are 
estimated from eOVi2. 

By the same procedure taken in Table 1 and Fig. 6, Table 
2 and Fig. 11 are prepared. Estimated values for BD2 are selected 
in the vicinity of point h' in Fig. 11. That is, eBC2 = -0.08, 
eBW2 = 1.11, ePH2 =1.06. In this case, the true value for BD2 
indicated as 'actual' is not within the hatched triangle. This 
is due to the poorly estimated parameter values for eBC3. However, 
an iterative operation described later can improve band 
parameter values for eBC3 . 
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Table 2 



Set of estimated value! x coordinate of cP22 e BW 2 ePH2A ePH2B 


EstBD2a 


0.360 


1.04 


o.yoz 


U.OOD 


EstBD2b 


0.365 


1.05 


1.02 


0.896 


EstBD2c 


0.370 


1.06 


1.05 


0.928 


EstBD2d 


0.375 


1.07 


1.09 


0.96 


EstBD2e 


0.380 


1.08 


1.12 


0.992 


EstBD2f 


0.385 


1.09 


1.16 


1.03 


EstBD2g 


0.390 


1.11 


1.2 


1.06 


EstBD2h 


0.395 


1.12 


1.23 


1.09 


EstBD2i 


0.400 


1.13 


1.27 


1.13 


EstBD2j 


0.405 


1.14 


1.31 


1.16 


EstBD2k 


0.410 


1.15 


1.35 


1.2 


EstBD21 


0.415 


1.17 


1.39 


1.24 


EstBD2m 


0.42C 


1.18 


1.43 


1.27 


EstBD2n 


0.425 


1.1S 


1.47 


1.31 
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(Step3) 

eOVi3 is found by band stripping eBD2 from 0V3 . D1-D2 plots 
regarding eOVi3 are shown in Fig. 12. As is clear from Fig. 12, 
since two local minima appear, it is found that there are at 
least three component bands within OV3 and the band center 
position could be estimated. Here, local minima Mini and Min 
3 indicate the existence of corresponding component bands BDi 
and BD3, respectively. 

When the positions of Max 1' and Max 3' corresponding to 
Mini (Xmini=-0 . 8257) and Min 3 (Xmin3=0 . 427 6) , respectively, are 
calculated, Xmaxi=-0 . 7757 and Xmax3=0 . 4403 . Min 3 and Max 3' found 
here are closer to the true value T3' compared to Min 3 and Max 
3' in Fig. 4 or 5. This is due to band stripping effect, which 
has reduced overlapping effect between component bands. 

Therefore, eBC3 in Step 1 is replaced by Xmax3 just obtained, 
and a set of eBD3 and its constraint conditions are recalculated 
in the same manner as in Step 1. At the same time, constraint 
conditions are also calculated from Fig. 12. The set of eBD3 
obtained by this and the constraint conditions are shown in Table 
3 and Fig. 13. According to the results, constraint conditions 
subjected to them are the narrower permissible regions. The 
intermediate point of points i' and j ' on SelLN2 is selected 
as a new estimated point . That is, eBC3 = 0.44, eBW3 = 0.81, and 
ePH3 =0.78. 
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Table3 



Set of estimated value 


X coordinate of eP32 


eBW 3 


ePH3A 


ePH3B 


r!iStJt>iJoa 




A CIO 


n /LDQ 




HiStx>UoD 


n 71 


u.DOD 


ft /ift 


ft 39fi 




ft 79 






0 372 


EstBD3d 


0.73 


0.683 


0.574 


0.422 


EstBD3e 


0.74 


0.706 


0.636 


0.476 


EstBD3f 


0.75 


0.730 


0.702 


0.535 


EstBD3g 


0.76 


0.753 


0.772 


0.598 


EstBD3h 


0.77 


0.777 


0.846 


0.666 


EstBD3i 


0.78 


0.801 


0.923 


0.739 


EstBD3j 


0.79 


0.824 


1.0 


0.816 


EstBD3k 


0.80 


0.848 


1.09 


0.899 


EstBD31 


0.81 


0.871 


1.17 


0.986 


EstBDSm 


0.82 


0.895 


1.26 


1.08 


EstBD3n 


0.83 


0.918 


1.36 


1.18 
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(Step4) 

When the band parameter values for BD2 are estimated in 
the same procedure as in the Step 2, eBC2 = -0.03, eBW2 = 1.04, 
and ePH2 = 1.01. 

When compared with the result of the Step 2, the suitability 
is improved. 

(StepS) 

Further, by iterating the same procedure as in the Step3, 
a new set of estimated values for BD3 is obtained. That is, eBC3 
= 0.449, eBW3 = 0.814, and ePHs = 0.800. 

(Step6) 

Further, band parameter values for BD2 are estimated and 
eBC2=-0.02, eBW2=1.04, and ePH2=0.967. Judging that estimated 
values are no longer improved with the iterative operation, the 
operation is hereby suspended. 
(Step7) 

Finally, band parameter values for BDi are estimated. 
This can be found by the geometry illustrated in Fig. 9. That 
is, eBCi = -0.642, eBWi = 1.15, and ePHi- 0.512. In addition, 
when the suitability of estimated parameter values for BD2 and 
BD3 is good, by D1-D2 plots where eBDi =OV3- eBD2 - eBD3, the band 
parameter values for eBDi can also be estimated. In Fig. 14, 
D1-D2 plots of eBDi and BDi are shown, whose shape is asymmetric 
and the curved portion of the left-hand side is smaller. 

These results show that estimated values of eBD2 or eBD3 
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are not good. Although the existence of the fourth component 
band cannot be denied, it cannot be judged from the estimated 
values obtained so far. Therefore, complementary matching is 
conducted assuming that there are three component bands. 

For confirmation, the estimated values at this stage are 
as follows: 

Estimated component band eBDi: eBCi = -0.642 

eBWi = 1.15 
ePHi = 0.512 



Estimated component band eBD2 : eBC2 = 0.002 

eBW2 = 1.01 
ePH2 = 0.967 

Estimated component band eBD3 : eBC3 = 0.4 49 

eBW3 = 0.814 
ePH3 = 0.800 



(Complementary Matching) 

The suitability of the component bands estimated at the 
present stage is examined. In Fig. 15, eBDi eBD2 eBD3 and eOV3 
are shown together with OV3. In Fig. 16, D1-D2 plots of eOV3 and 
OV3 are shown. In the figures, the iso-wavenumber lines between 
OV3 and eOV3 were also drawn and the profiles of eOV3 and OV3 are 
not coincident. However, since it is hard to examine from Fig. 15 
or 16 which component band has good or poor suitability for an 
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estimation, the improvement of band parameter values should be 
performed and evaluated according to other criteria. 

Therefore, with a complementary estimated component band, 
the improvement of band parameter, values is achieved. Here, 
a complementary estimated component band is the component band 
subtracting all the other component bands except only one 
component band between estimated component bands from the 
spectrum. That is, the complementary estimated component band 
cBD3 with respect to BD3 is represented as cBD3 =OV3-eBDi-eBD2. 
If estimated values of all the band parameters are coincident 
with the true values, eBD3 and cBD3 are also coincident . Likewise, 
the complementary estimated component band CBD2 or cBDi can be 
calculated. 

Figure 17 shows D1-D2 plots of eBD3 and cBD3. Likewise, 
Figs. 18 and 19 show D1-D2 plots of eBD2 and cBD2, and eBDi and 
cBDi, respectively. In Fig. 17, since Min3 (Xmin3=0 . 4363 ) , the 
characteristic point of cBD3, is close to Max3' (Xmax3=0 . 44 65) , 
the band center position may be estimated at the corresponding 
local maximum or minimum value of eBD3. The right-hand side of 
the plot is smaller compared with eBD3. This may be because eBD2 
or eBD3 is estimated smaller than the true value. 

Next, in Fig. 18, eBD2 is examined. Since Min2 (Xmin2=- 
0.0195), the characteristic point of cBD2, is close to Max2' 
(Xmax2=-0. 0152) , the band center position of BD2 may be estimated 
at the corresponding local maximum or minimum value of eBD2. 
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The plot of a single band is heart-shaped since the peak height 
ePH2 is estimated smaller than the true value. 

Further, in Fig. 19, eBDi is examined. Mini (Xmini=- 
0.6390) , the characteristic point of cBDi, and Maxl' (Xmaxi=- 
0.7038) are not close to each other. In addition, since the 
degree of the symmetry of the shape of cBDi is not good, the 
estimated values for the band parameters are not good. 

Now that suitability of an estimation of each component 
band was made clear, suitability is enhanced by adjusting each 
band parameter value. As the criteria for suitability of an 
evaluation value, SumLS, the total sum of the distances of the 
iso-wavenumber lines between OV3 and eOV3 on D1-D2 plots is 
provided and by adjusting each parameter value so as to make 
the evaluation value small, suitability is enhanced. 

As examined from the profiles of Figs. 17 to 19, first, 
parameter values for eBDiare adjusted. Adjustment is performed 
in the order of eBDi, eBWi, and ePHi. New eBCi is determined as 
the estimated value when SumLS is the minimum by changing it 
in the vicinity of the current eBCi. When SumLS is calculated 
and eBCi is adjusted regarding the total range of X (-2 < X <2) , 
eBCi=-0.576. When SumLS is calculated and eBCi is adjusted 
within the limited range of the inflection points P11 to P12, 
eBCi=-0.594. Further, when SumLS is calculated and eBCi is 
adjusted within the limited range of the secondary inflection 
points Q12 to Q13, eBCi— 0.600. Since the true value of BCi is 



40 



-0.5978, at the early stage of the optimization of band center 
position, that is, at the stage where suitability of the 
estimation is poor, it is better to estimate the band center 
position limiting to the vicinity of a corresponding local 
maximum or minimum to calculate SumLS and adjust the parameters. 

Likewise, as for estimating eBWi, adjustment is performed. 
As for estimating ePHi, it is better to find such ePHi as making 
SumLS minimum regarding the entire ranges of X. 

The same kind of operation is conducted on eBD2 and eBD3 
as well. Figures 20 to 22 show D1-D2 plots of the complementary 
estimated component band after iterating a series of the 
parameter values adjustment three times. Compared Figs. 17 to 
19, suitability of the band parameter values is improved. 

The optimization of all the band parameter values cannot 
be simultaneously achieved when the estimated value of each band 
parameter values is adjusted alone. The optimization of the 
estimated values for eBDi are so adjusted as to enhance the degree 
of coincidence between eBDi and cBDi as well as to improve the 
symmetry of cBD2 which is the adjacent complementary estimated 
component band in Fig. 18. 

Here, the criteria of symmetry are so defined that two 
inflection points are symmetrical to the Y-axis. That is, for 
example, the parameter values are so adjusted that the length 
of the line segment OPu is equal to that of the line segment 
OP12. 
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For example, in adjusting ePHi, ePHii where SumLS is the 
minimum (ePHii = 0. 64 6) is found. Next, ePHij is found where the 
length of line segment OP11 is equal to that of line segment OP12. 
Then its mean value (ePHii+ ePHij)/2 is selected as an optimum 
estimated value. The reason for selecting the mean value is that 
the symmetry of CBD2 is affected from the band of both sides. 
By such iterated operations, ePHii and ePHij become coincident. 

In adjusting the band center position eBCi, the above 
operation is not always required since the good estimated value 
is obtained for eBCi where SumLS has the minimum value. Likewise, 
the improvement of the coincidence between eBD2 and CBD2 and in 
the symmetry of CBD2 is achieved. The results are shown in Figs. 
23 to 25. Finally, the good results were obtained as follows: 
Estimated component band eBDi: eBCi = -0.597 

eBWi = 1.094 
ePHi = 0.648 

Estimated component band eBD2 : eBC2 = 0.004 

eBW2 = 0. 952 
ePH2 = 1.070 

Estimated component band eBD3 : eBC3 = 0.456 

eBW3 = 0.794 
ePH3 = 0.800 

Finally, the examination of the existence of the hidden 
component bands or the evaluation of the residual component rBD 
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is performed, rBD is represented as rBD = OV3-eBDi-eBD2-eBD3. 
Virtual component band vBD whose parameter values are already 
known is introduced and D1-D2 plots of the data with rBD, 
remaining component bands, and vBD added are evaluated. Since 
the derivative plot of vBD is already known, it can be evaluated 
by the shape of the overlapping residual component. In 
addition, an evaluation of a complementary estimated component 
band is also made. Thus, when a hidden component band can be 
estimated, it is a good way to return to the first step, finding 
the estimated value of parameter values for a new component band 
from constraint conditions, thereby improving the parameter 
values of the said component band by the similar method. 

The estimation method of the present invention was 
described above by using the simulation data. In estimating 
actually measured data, it is necessary to differentiate 
digitally in order to obtain derivatives plots. Therefore, in 
the case where the band to be measured is sharp, it is necessary 
to measure with the small data interval for measurement. 

Particularly, in the present invention, as for the 
analyzed object, what matters is a method of analyzing the 
obtained data and in order to obtain the data, the conventional 
spectrophotometer can be used as it is thereby requiring no new 
spectrophotometer. The component bands can be estimated from 
the spectral data of the analyzed object by outputting spectral 
information, thereafter putting into a data station or a computer 
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with programs developed based on the proposed method. These 
programs can be installed into the spectrophotometer. 
(Effect of the Invention) 

When the procedure of the spectral analysis method of the 
present invention is used, several component bands can be easily 
estimated as for the analyzed object having a spectral profile 
which contains several component bands. 

Summary on the procedure of the spectral analysis 
described above is as follows: 

(1) The band center position of corresponding component 
band is determined noting the typical local minimum, and/or local 
maximum of the two-dimensional derivative plot (including D1-D2 
plots) of the spectra of the analyzed object. 

(2) The candidates for the bandwidth values are determined 
by selecting candidates of several points in the vicinity of 
zero-crossing points of the second derivatives, estimating the 
bandwidth of the corresponding component band from the 
candidates of inflection points, and estimating the peak height 
of the corresponding component band from the algebraic 
geometrical condition of two-dimensional derivative plot. 
Further, band parameter values are imposed by the constraint 
conditions and the estimated value is determined. 

(3) The band parameter values are estimated by finding 
the component bands in order by band stripping. 

(4) The estimated band parameter values are optimized and 
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improved by using the complementary matching method and 
considering the overlapping effects among component bands, 

(5) The band parameter values for component bands are 
determined among best fitting band shape since the profile of 
each component band must be line symmetry when each component 
band is extracted. 

(6) Suitability of the band shape and existence of hidden 
bands are examined by observing profile of each component band 
after optimization of the band parameter values. 

Description of the preferred embodiment described herein 
is illustrative and not restrictive, the scope of the invention 
being indicated by the appended claims and all variations which 
come within the meaning of the claims are intended to be embraced 
therein . 
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